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Abstract 

We generalize the Weisskopf-Wigner theory for the hne shape and transition rates of decay- 
ing states to the case of the energy-driven stochastic Schrodinger equation that has been used 
as a phenomenology for state vector reduction. Within the standard approximations used in 
the Weisskopf-Wigner analysis, and assuming that the perturbing potential inducing the decay 
has vanishing matrix elements within the degenerate manifold containing the decaying state, the 
stochastic Schrodinger equation linearizes. Solving the linearized equations, we find no change from 
the standard analysis in the line shape or the transition rate per unit time. The only effect of the 
stochastic terms is to alter the early time transient behavior of the decay, in a way that eliminates 
the quantum Zeno effect. We apply our results to estimate experimental bounds on the parameter 
governing the stochastic effects. In an Added Note, elegant stochastic-theoretic methods suggested 
by Diosi are used to rederive the principal results, without the assumptions needed to linearize the 
stochastic equation, and to give analogous results for the Rabi oscillations of a two-level system. 
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I. INTRODUCTION 



There has recently been considerable interest in the possibility that quantum mechan- 
ics, and the Schrodinger equation, may be modified at a very low level by effects arising 
from Planck scale physics. Such speculations have been motivated on the one hand by 
considerations from string theory and quantum gravity [Q, and on the other hand by 
efforts 0, 1^, Q to achieve an objective equation describing state vector reduction. The 
majority of the objective reduction discussions fall into two classes: those that postulate 
a stochastic process producing spatial localization [^, and those that postulate an analo- 
gous stochastic process leading to localization in energy |Q (the so-called "energy-driven" 
stochastic Schrodinger equation.) Both the spatial localization and the energy localization 
stochastic Schrodinger equations avoid problems with superluminal signal propagation that 
characterize attempts at deterministic nonlinear modifications of the Schrodinger equation 
1^. We find the energy-driven approach particularly appealing because it is energy conserv- 
ing, leads with no approximations to Born rule probabilities and to the Liiders projection 
postulate, has sensible clustering properties, and when environmental interactions are taken 
into account explains state vector reduction with a single Planck scale stochastic parameter 

Although physical prejudices might suggest a Planck scale magnitude for the stochastic 
parameter in the energy-driven equation, one can instead take the point of view that the 
stochastic parameter can have a priori any value, and use current experimental information 
to place bounds on it. This approach has been pursued [0 in the context of particle physics 
systems that exhibit oscillations between different mass eigenstates (the i^'-meson, 5-meson, 
and neutrino systems), with results that are summarized in the final section of this paper. 
An alternative source of bounds on the stochastic parameter could come from experiments 
observing decays and line shapes in atomic and particle systems, since if the stochastic 
terms in the Schrodinger equation were to change the standard Weisskopf-Wigner analysis 
of decay processes in a significant way, then observable effects could result. Thus, to pursue 
phenomenological studies of the energy-driven equation, it is important to generalize the 
standard Weisskopf-Wigner decay theory P] to include effects of the energy-driven stochastic 
terms. This is the problem that is analyzed in this paper. 
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II. THE ENERGY-DRIVEN STOCHASTIC SCHRODINGER EQUATION AND 
PROPERTIES OF THE ITO STOCHASTIC CALCULUS 



Letting denote a unit normalized Schrodinger picture state vector, the standard form 
0, ^, P of the energy-driven stochastic Schrodinger equation is (with h= 1) 

d\ij) = -iH\ij)dt-\a\H-{H)f\ij)dt + \a{H-{H))\ij)dWt . (la) 
o 2 

Here H is the Hamiltonian, (H) = {iIj\H\iIj) is the expectation of the Hamiltonian in the 

state lip), 0" is a numerical parameter governing the strength of the stochasticity, and dWt 

is an Ito stochastic differential that, together with dt, obeys the standard Ito calculus rules 

i 

dW^ = dt , dWtdt = dt'^ = . (lb) 

By construction, the nonlinear evolution of Eq. (|T^ guarantees the preservation in time of 
the unit normalization of the state vector 

In the following sections, we shall need a number of properties of the Ito calculus that 
we summarize here. First of all, in the Ito calculus the Leibnitz chain rule generalizes to 

d{AB) = {A + dA) [B + dB) - AB = {dA)B + AdB + dAdB , (2a) 

with the final term dAdB contributing a term proportional to dt when the dWt contributions 
to both dA and dB are nonzero. Letting Wt be the Brownian motion 

Wt= I dWu , (2b) 
Jo 

we see in particular that 

dex.p{aWt) = exp{aWt)[exp{adWt) - 1] = exp{aWt)[adWt + ^a'^dt] . (2c) 

Letting E[...] denote the stochastic expectation of its argument, and letting A(t) denote any 
function of the stochastic process up to time t, we have 

E[dWtA{t)] = , (3a) 

since the Ito differential refers to the time interval from t to t + dt, and hence is statistically 
independent of the process up to time t. Thus, taking the expectation of Eq. (^, we get 
the differential equation 

dE[exp{aWt)] = E[exp{aWt)]^a^dt , (3b) 
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which can be immediately integrated to give 

E[exp{aWt)] = exp{^aH) , (3c) 

a result that will be needed later on. 

Let us make an elementary application of the Ito formalism, to write the stochastic 
Schrodinger equation of Eq. (p!a]) in an equivalent form. First of all, forming the density 
matrix 

P=l^)(^l , (4a) 

we have from Eq. (pa|), 

dp={d\^)){ij\ + \ij)d{ij\+d\ij)d{ij\ , (4b) 



which on substitution of Eq. ( pTa]) and use of the Ito calculus rules of Eq. ([Tbl) gives the 
evolution equation for the density matrix, 

dp = t[p,H]dt-^a'[H,[H,p]]dt + ^a[p,[p,H]]dWt . (4c) 

Taking the stochastic expectation of this equation, using Eq. (Pa]), gives a differential equa- 
tion of the Lindblad type 0] for E[p], 

^ = ^[E[p],H]-la'[H,[H,E[p]]] . (5) 

The fact that this equation is linear (in contrast to Eq. (^, which is nonlinear) is the 
fundamental reason |^ why Eq. ( p!aD does not give rise to superluminal signal propagation. 



III. INITIAL FORMULATION OF THE DECAY PROBLEM 

Let us now formulate the decay problem for the stochastic Schrodinger equation of 
Eq. (p!aD, following the standard procedure for the usual Schrodinger equation without 
stochastic terms. We suppose that for times t < the Hamiltonian H is given by an 
unperturbed Hamiltonian Ho, with eigenstates \n) and eigenvalues En, 

Ho\n) = En\n) , (6a) 

and that the system under consideration is in an eigenstate \sa) with eigenvalue Eg, which 
is one of a set of degenerate energy eigenstates \sa) , a = 1,...,D. Because Eq. ( p!a]) 
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acts as an ordinary Schrodinger evolution on a state IV') that is an energy eigenstate, the 
system remains in the state \sa) as long as the Hamiltonian remains equal to Hq. Hence the 
starting point for the standard decay analysis |11| is also a consistent starting point for its 
stochastic extension under Eq. (|T^). As in the standard procedure, we assume that at t = 
a time-independent perturbation V is switched on, so that for times t > the Hamiltonian 
is H = Hq + V . The initial state \sa) is then no longer an energy eigenstate, and so will 
decay into various other states |m); our problem, as in the usual case, is to find the partial 
transition rates for this decay and the probability amplitude for the system to remain in the 
initial degenerate group of states. 

In formulating this problem, it is convenient to expand the state {ip) over the basis |n) 
and, as in the standard case, to remove the Schrodinger time evolution associated with the 
unperturbed Hamiltonian Hq, by writing 



Substituting Eq. (|6bD into Eq. (p!a]) , and projecting on (m|, it is a matter of straightforward 
but somewhat tedious algebra to compute the stochastic evolution equation for Cm(t), with 
the result 



n 

am = -IcT^Em - {H)fdt + \a{E^ - {H))dWt , 
o Z 

Pmn = -iVmneXp[i{Era - En)t]dt - \a\Em + En - 2{H))Vmn + {V'^)mn\ exp[i{Em - E„ 

o 

+ -(rVrnneyip[iiEm- En)t\dWt . 
The corresponding expression for (H) is 

{H) = Y,En\Cn{t)\' + J2VrnneXp[l{Ern-En)t]C*Jt)Cn{t) . (7b) 
n mn 

In these equations Vmn and (V^^)mn denote the respective matrix elements 

Vmn = {m\V\n) , {V^)n.n = {m\V^\n) . (7c) 



IV. APPROXIMATION TO LEADING ORDER IN V 



Equations d^-c) are a complicated, nonlinear set of stochastic differential equations, and 
so to solve them approximations will be needed. Following the Weisskopf-Wigner analysis, 



we shall make the approximation of regarding as a small perturbation. The coefficients 
Cm, 1TI ^ {sa} for states not in the initial degenerate manifold will then be of order 0{V), 
and we neglect 0(V^^) and higher contributions to them (except those arising implicitly 
through our solution for the Cg^)- On the other hand, the coefficients Cs^ of states in the 
degenerate manifold can be of order unity, and we calculate these coefficients to order 
accuracy, neglecting corrections of order and higher. In a similar fashion, in expressions 
involving the stochasticity parameter a, we shall retain terms of order aV and its powers 
(aV)"^, etc., but shall neglect terms of order aV"^ and higher that involve extra factors of V 
relative to the terms that we are retaining. Finally, although we shall see that Em — Eg is 
effectively small, we shall retain all terms of order a{Em — Eg), a'^{Em — Eg), [a{Em — Eg)]"^, 
etc., but shall drop terms (J {Em — Eg)aO{V'^) that are smaller than these by a factor of 
order aV"^ or V"^. 

Making use of the perturbative ordering of the coefficients C„, we begin by simplifying 
and approximating the expression in Eq. ( [TED for {H). Separating off the states in the initial 
degenerate manifold, the sum in Eq. ( [TED becomes 

{H) = EgY,\Csf + Y.^s.^^C:Cg, + 0{y^) . (8a) 

a ab 

However, since the state vector IV') remains unit normalized, we have 

El^-I' = l- E |C™P = 1 + 0(^2) , (8b) 

and so we have 

{H) = Eg + Y,Vs.s,ClCg, + 0{y') . (8c) 

ab 

If we substitute Eq. (|c|) back into Eq. (|7aD, we are still left with a nonlinear set of 
equations. Therefore we shall also introduce the simplifying assumption that the perturbing 
potential has vanishing matrix elements within the degenerate manifold containing the initial 
state, so that 

Vg^g, = 0, a,6 = l,...,D . (9a) 

There are important, physically relevant cases in which Eq. (^) is obeyed as a result of 
selection rules. For example, for radiative decays treated in the electric dipole approximation, 
with Ho taken as the atomic Hamiltonian plus the free radiation Hamiltonian, and with V 
taken as the atomic couphng to the transverse electromagnetic modes, Eq. ( pa]) is obeyed as 
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a result of parity invariance when the states in the initial degenerate manifold all have the 
same parity. (We caution, however, that Eq. ( pa]) is not valid for the analysis of K (or B) 
meson systems when Hq is taken as the strong interaction Hamiltonian. Thus here either 
one has to employ the nonlinear equations following from Eq. (Pcj), or one has to redefine 
Hq so as to impose Eq. (^) by including in Hq the \^S\ = 2 (or |AC| = 2) weak interaction 
effective Hamiltonian terms, with V defined to contain only the | AS"! = 1 (or | AC| = 1) weak 
interaction terms responsible for K (or B) meson decays. Such a redefinition is consistent 
in the vacuum saturation approximation for the \AS\ = 2 (or |AC| = 2) terms.) 
With the simplifying assumption of Eq. (|9a|) , Eq. (|c|) becomes simply 

{H) = E, + 0{V') . (9b) 

Substituting this into Eq. (^), and dropping terms that are not of leading order in V in the 
sense defined above, Eqs. (|7alJ7BD simplify to the following set of linear equations, 

n 

l^^n = \(yy^n, l^± = -lVrnn-\A{Ern + En-2Es)V^n + {V^)rnn] • (10a) 

Corresponding to the magnitude ordering of the coefficients Cm introduced above, it is 
convenient to rewrite Eq. (|10a]) as separate equations for the two cases, m G {s^} and 
m ^ {Sa}- For m G {Sa} the coefficients al ' vanish; separating the sum over n into terms 
where n G {s^} and n ^ {sa}, using the assumption of Eq. (|9al) , and dropping terms of 
nonleading order in V, we get 

^ b 

+ J2 ^MKEs - Er^m^ilidWt + 7i%dt)Cn{t) , 
7£l = l^Vs^n , li% ^ -iVsanfn , (10b) 

where we have introduced the definition 

fn = l-lcj\En-Es) . (10c) 



7 



For m ^ {sa} the coefficients am are nonzero, but only the terms with n G {sa} have to 
be retained in the sum over n, and so we similarly get 

a 

7^1 = IcrVms. , 7^1 ^ -^VmsJm ■ (lOd) 

Equations (|IO|a-d) are the basic system of stochastic differential equations that we shall solve 
in the subsequent sections. 



V. EQUATIONS FOR EXPECTATIONS OF THE COEFFICIENTS 

The principal quantities that we wish to calculate are the expectations -E[|Cm(t)P] of the 
squared magnitudes of the coefficients, since these give the expectations of the probabilities 
for the various states to be occupied. We shall show in this section that, within our ap- 
proximations, these can be directly related to the expectations E[Cm(t)] of the coefficients 
themselves, for which we shall derive a closed, linear set of ordinary differential equations. 

Again, we consider separately the cases m G {sa} and m ^ {sa}- For Cg^, we write 

C,St) = E[C,M + ^a{t) , (11a) 

with E[Aa{t)] = 0, and with Aa{0) = since the stochastic terms in the differential equation 
act only after t = 0. However, referring to Eq. ( |10b| ) we see that dCs^ is of order V"^, and so 
Aa{t) must also be of order V^. Therefore 

E[\CsAm = \E[CgM\' + 0{V') , (lib) 

and so to the accuracy to which we are working, we can compute i^^[|C5^()f:)p] from the 
expectation E[CsSt)], ignoring the effects of fluctuations. 

We consider next i?[|Cm(t)p] for m ^ {Sa}- Applying the Ito rule of Eq. (|2a|) , we have 

dE[\CUt)\'] = E[{dC*^{t))Cm{t) + C*„MdCUt) + dC*^{t)dCUt)] ■ (12a) 



Substituting Eq. ( |10d|) for dCmif) and using Eq. (|3a|) , which eliminates the dWt contribu- 



tions, and using the fact that to leading order in V we can replace Cg^ {t) by its expectation, 
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we get after some algebraic simplification the formula 

a 

- exp[z(E„ - E,)t]tr^E[C:M E VmsMCsSt)] 

a 

+ \cr'\J2^^msMCsSm' , (12b) 

a 

which can be integrated to give once the expectations E[Cm{t)] and E[Cs^(t)] 

are known. 

To get a closed set of equations for the expectations of the coefficients, we simply take 
the expectations of Eqs. ( |lUb| ) and ( p^Odl ), and use Eq. (^), which again eliminates the dWt 



contributions. For -^[Cs^] we thus get 



dt 

b 



eMKEs - En)t\{-i)Vs^nfnE[Cn{t)] , (13a) 

n^{sa} 



while for E[Cm{t)\ with m ^ {sa} we find 
d . 1 



j^E[C^{t)] = -^a\Era-E,f)E[C^{t)] 

+ eM^{E^- Es)t]Y{-i)VmsJmE\CsM ■ (13b) 

b 

VI. SOLUTIONS FOR EXPECTATIONS OF THE COEFFICIENTS 

We proceed now to solve the linear system of equations for the expectations of the coeffi- 
cients given in Eqs. (|13|a,b). Since the problem is defined on the half line t > 0, the natural 
way to do this is by using the Laplace transform. Defining 

/"OO 

gUp) = / dtexp{-pt)E[CUt)] , (14a) 
Jo 

we have, by an integration by parts, 

rftexp(-pt)^^i^M = pg^^^p) _ E[CUO)] , (14b) 



and also 



POO 

/ dtexp{-pt)exp[i{Em- E„)t]E[Cnit)]= gnip-iEra + iEn) , (14c) 
Jo 



with the integrals in Eqs. (0a-c) defining analytic functions of p in the right hand half plane 
Rep > 0. The inversion of the Laplace transform is given by the formula 

E[CUt)] = TT- / dpexp{pt)gmip) , (14d) 

with e > an infinitesimal positive constant. 

Taking the Laplace transform of Eqs. (|13|a,b), and using the initial conditions E[Cs^{0)] = 
CsM = SaA and E[C„,{0)] = C„(0) = 0,m ^ {s,}, we get 

b 

+ ^ {-i)Vs^nfn9n{p + 'i'En-iEs) , (15a) 

n^{sa} 

and for m ^ {s^}, 

P9m{p) = -^CX^(^m - EsY)gm{p) 

+ y^X-iWmsJmgsjiP + iEs - lEm) ■ (15b) 

Solving Eq. (|15bD for gm{p), Tn ^ {sa\, and shifting p — > p + iEm in the solution, we get 



gm{P + iEm) = [P + + ^0-^(^m - Esf] ^ ^{-iWmsJmQsSP + iEs) ■ (16a) 



b 



Shifting p ^ p + iEs in Eq- ( |15a| ), and then substituting Eq. ( |16a| ), we get an algebraic 



equation for the set of quantities gstip + iE. 
^Kabgsiip + iEs) -- 



aA 1 

2 



K.b = {p + ^Es)6^, + la\V')s^s^+ ^ ■ /"f v^^"^" ^ w • (16b) 

8 7f ,P + ^Es + l{Em - Es)fm 

In physically interpreting these equations, we must remember that the Laplace transform 
variable p is related to the usual energy variable E hj p = —iE. Making this substitution 
in Eqs. (|T6|a,b) we have respectively 

gm{-iE + iEm) = [-iE + lEm + ^(J^{Em - ^.)T' Y^{-i)VmsJ^gs,{-iE + iEs) , (17a) 

b 

and 

Kabgst{-iE + iEs) = SaA , 

b 

Kab = {-tE + tEs)6ab + -a\V^)s s,+ y —p, ^i^^^''";!r^'' .(17b) 



-iE + tEs + i{Em-Es)f„ 



10 



Corresponding to the changes of variable that have been made, the inversion formulas become 
E[CsSt)] = ^exp(zE,t) / dEexp{-iEt)gsS-iE + tEs) , 

^'^ Jie-OQ 

E[C„,<^{sa}{t)] = —exp{iE^t) dEexp{-iEt)g^{-iE + iEm) ■ (17c) 

Inspecting the equation for the kernel Kab, we see that apart from order V'^ terms it 
is a diagonal matrix {—iE + iEs)6ab- Hence the solution gs^{—iE + iEs), on the inversion 
contour of integration, will be appreciable only in the vicinity of E = ie + Eg, that is, only 
near energy shell. This motivates the Weisskopf-Wigner approximation of replacing E in 
the denominator of the final term in Kab by ie + Es, with the result that Kab then becomes 
a linear function of E. Before making this approximation, the kernel Kab has a non-trivial 
dependence on the stochasticity parameter a. However, after making the Weisskopf-Wigner 
approximation, this cr dependence completely cancels: 

rrv^{sa} 

1__2/ta2'\ I \ ^ fm^SamVmSf, 



-lE + iEs + i{Em- Es)f„ 

m^{Sa} 



1 2n/2\ I \^ [1 - {i/9>)a'^{Em - Es)\Vs^rnVmsi, 

S"" '-^'^^ e + ^{E^-Eg) 

E^, - E. - ' ^^^^^ 



m^{Sa} 

where in the final step we have made use of the condition of Eq. (Pa|). Thus in the Weisskopf- 
Wigner approximation, the kernel Kab appearing in Eq. ( |17b|) simplifies to 

Kab = HE + tEg)Sab + Yl "'^""^^"^^^ 



-, Em — Eg — ie 



-iE + iEg)5ab + iMab + ]^Tab 



Mab = 2^ ^~E^ 



1 Eg — Em 



Tab = 27r 5^ Vg^mVmsAEm ' Eg) , (I8b) 
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with P in the definition of the "mass matrix" Mab the principal value. These are the same as 
the formulas for the kernel in the absence of the stochastic terms in the Schrodinger equation. 
Thus, in the Weisskopf-Wigner approximation, the solution for E[Cs^(t)] is unmodified by 
the stochastic effects, and hence the Lorentzian line profile and the decay rate of the state 
are unaffected by the a terms. 

The solution for E[Cm(t)] with m ^ {sa} does retain a dependence on the stochastic 
parameter. To study this, let us specialize to the case D = 1 of a non-degenerate initial 
state. The expression in Eq. (|18b| ) for the kernel now becomes the 1x1 matrix 

K{E) = -iE + iEs + iM+^T , (19a) 

with M and F real numbers given by 

M = \^ p VsmVms 

^ E - E 

V = 2TlY,ysny^s5{Em-E,) . (19b) 

Thus, Eq. ( p.7b|) has the immediate solution 

gs{-iE + iE,) = K{E)-^ , (19c) 

which when substituted into Eq. (|17aD yields 

gm{-iE + lE^) = [~iE + lE^^ + \(j\Em - Esf]-W^sfm{E -Es-M+ '-T)-^ . (19d) 

o Z 

Substituting these equations into the inversion formulas of Eqs. (|17c|) , and doing elemen- 
tary contour integrations, we find 

E[Cs{t)] = exp{-iMt - ^rt), (20a) 
E[C„.Mt)] = _ j^^l'^ _ (^^M^iE„, -Es- M)t - ^Tt] - exp[-^a'iE„, - E,)H]^ . 

/^From Eq. (|llb| ) we thus get 

E[|C,(t)|2] =exp(-n) , (20b) 



which identifies T as the transition rate per unit time out of the initial state. Finally, 
substituting Eq. ( |20a|) into Eq. ( |12b| ), simplifying to leading order in V, and integrating 
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with respect to t, we get 



- 2 exp[-ia2(E„ - E,)H - ht] cos[(^, ~E^ + M)t]^ . (20c) 

This completes our solution for the expectations of the coefficients, and their squared magni- 
tudes, in the case of a non-degenerate initial state. We see that after a time t large compared 
with the lifetime T~^, we obtain 

exhibiting the standard Lorentzian profile with no dependence on the stochasticity parameter 
a. 



VII. SMALL TIME AND GOLDEN RULE APPROXIMATIONS 

Let us now study the behavior of Eq. (|20c|) for small and large values of the time t. Since 
within our approximations we have a'^{Em — ^ o''^[{Es — E^ + M)^ + |r^], we can 
rewrite Eq. ( P0c| ) as 



|2 



-2exp[-^a^{{Es - Em + Mf + ^V'')t - ^H] cos[(E, -Em + M)t]^ . 

In the limit as t ^ 0, we can develop the exponential and cosine functions in power series 
expansions, giving 

E[\Cm^M'] ^ iKnsP Q^'t + + 0{ea\E, - Emf) + 0{e)^ . (21b) 

Thus the leading small time behavior of the summed expected probability in the decay 
channels is 

J2\Vms\'\^'t={V'),,^aH , (21c) 

where in evaluating the sum we have employed the condition of Eq. (|9^). We shall verify 
this result by another method in Sec. |X[ where we discuss its implications for the quantum 
Zeno effect, and in Sec. ^ shall apply it to estimating bounds on a. 
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Let us next consider the large time behavior imphed by Eq. (|21a|) . Once t is large 
enough so that \{Em — Es)t\ is large for all energies Em not infinitesimally close to Eg, 
we can evaluate the summed expected probability in the decay channels by making the 
"golden rule" approximation [|T^. This approximation treats the factors multiplying 
in Eq. ( [^la| ), which are sharply peaked around Em = Eg, as if they were equal to a Dirac 
delta function of strength given by the integral of these factors over energy. We then have 



- 2exp[-i(7'((E, -Em + Mf + h:^)t - h:t\ cos[(E, -Em + M)t]^ 

^ Y,\yms?KEm-Es) 

r . ^^ exp(-n) + 1 - 2exp[-|a^((AE)^ + \T^)t - in] cos[(AE)t] 

" i-oo'^^^^ (Kwrw 

= ^n^vm,t] , (22a) 

with r as given in Eq. (|19bD and with the function F[y4,t] defined by 



oo 



exp(-n) + 1-2 exp[-A{u^ + IVH^) - |n] cosm 



To evaluate we note that [|l3l 

^ ^ /-"^ exp(-rt) + 1 -2exp(-irt)cosM 27r, , , , 

F[0,t] = jju-^ ^,^1^,^, ' = -[l-exp(-n)] , (23a) 



and 13 



dF[A,t] 
dA 



f°° 11 

/ ciM2exp[-v4(u^ + -T^t^) - -Tt] cosm 

J-oo 4 2 

27i^/^exp{-^Tt-^ATH^)A-^/^exp[-l/{4A)] . (23b) 



Thus, integrating Eq. (|23b|) with respect to A we get 



F[At] = ^[1 - exp(-rt)] + C[A,t] , (23c) 



with the correction term C[A,)f:] given by 



C[A,t]=47T^/^exp{--Tt) [ dvexp[~-{v^TH^ + l/v^)] . (24a) 

2 Jo 4 
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Since the exponentials of negative arguments in Eq. ( p4a| ) are bounded by their maxima 
over the range of integration, we have 

\C[A,t]\ <A7r^/^A^/^exp[-l/{AA)]=2a{^y/^exp{-2t/a^) . (24b) 

So when Ft is of order unity, the contribution of the correction term C[y4, t] is of order 
caF^/^ exp[— 2/(cr^r)] ~ c'aV exTp[—c" / (aVY], with c,c',c" constants, which is exponen- 
tially small and can be neglected in our approximation scheme. Thus we are justified in 
approximating 

F[A,t] ^ F[0,t] = ^[1 - exp(-n)] , (24c) 
which when substituted back into Eq. (|22a] ) gives 

^E[|C„^,(t)|2] = l-exp(-n) = l-|C,(t)p , (24d) 

verifying that the approximations used in our calculation are consistent with maintenance 
of the unitarity sum rule (the unit state vector normalization condition). 

VIII. SOLUTION TO THE STOCHASTIC EQUATION FOR C^^, 

Since we see from Eqs. (M ) and (|OcD that E[|C„-,,(t)|2 differs from \E[C„^^s{t)]\'^ , the 
stochastic fluctuations in Cm^sit) are evidently playing a role. Let us now demonstrate this 
directly by solving the stochastic differential equation for Cm^s(t). Specializing to the case 
of a non-degenerate initial state, approximating Cs(t) ^ E[Cs(t)], and using Eq. ( P0a| ) for 
E[Cs{t)], Eq. (|T0g|) becomes 



+ exp[z(E^ - - M)t - ht] i^y;^JWt + itldt) , 

= \<y{En. - E.) , = -la^(E„ - E^f = -1(««)^ , 

7i'l = l^Vrns , 1^2 ^ -^Vrnsfm • (25a) 

For general values of the coefficients am and fm's >, Eq. (|25|) can be integrated by using 
Eqs. (0a-c) to find a stochastic integrating factor for the Cm terms (see the Appendix), with 
the result 

Cm{t) = exp[««W^i - (««)2t] j\xp[z{Em-E,-M)u-^Tu-ag^W^+{ag^fu] 

X b2dW^+{jl^l-agh2,)du] . (25b) 
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Using Eqs. (0a-c), it is easy to verify directly that Eq. ( |25b|) solves Eq. ( |25a|) . If we now 
examine Eq. ( p5b| ) more closely, using the specific expressions for the coefficients am'^^ and 
7ms given in Eq. ( p5a| ), we find that within the approximation of neglecting terms of relative 
order aV'^, the integrand in Eq. (|25b|) is an exact stochastic differential. Thus the integration 



can be carried out explicitly (see the Appendix), with the result 

Cm^sit) = — (25c) 

exp[^(Em -Es- M)t - ht] - exp[^or(Em - E,)Wt - ^cr\E^ - E,)H] 

This expression can be easily verified, by use of Eqs. (2a-c), to be the solution to Eq. ( |25a] ) 
(up to a residual error of relative order aV"^). Using Eq. (|3^) to take the expectation of 
Eq. ( p5c| ), we recover the result of Eq. (|20a|) . From Eq. ( p5c|) we find an explicit formula for 

W P 



|Cm^s(^)| 



2 



{E, -E^ + MY + 



exp(-n) + exp[a(Em - Es)Wt - ^a^E^ - E,)h] (25d) 
- 2 exp[^a(Em - E,)Wt - ^a\E^ - E.ft - ^Et] cos(Em -E^- M)t 



Again using Eq. ( p^ to take the expectation of this formula, we recover the result of 
Eq. O. 



IX. STOCHASTIC SUPPRESSION OF THE QUANTUM ZENO EFFECT 



In Eqs. (|2T[b,c) we saw that the leading small time behavior of the summed expected 
probability in the decay channels is 

iy')ss\aH , (26) 

rather than the result (y'^)sst^ that would hold for vanishing a. As a result, i?[|Cs(t)P] — 1 
vanishes linearly in t for nonzero a, rather than quadratically in t as for the unmodified 
Schrodinger equation. Since the quadratic vanishing of |Cs(t)p — 1 in standard quantum 
mechanics is the origin of the quantum Zeno effect |T^, we conclude that in the energy 



driven stochastic Schrodinger equation, the quantum Zeno effect is suppressed. 
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Let us verify this directly from the stochastic differential equation of Eq. (|T^), in analogy 
with the direct calculation [|15] of |Cs(t)P — 1 for small times for the ordinary Schrodinger 
equation. Applying the Ito rule of Eq. (|2a|), we have 

d\{smsm\=o = {smd\smt=o+{smd\sit)n=o+{smd^^ ■ 

(27a) 

From Eq. (^) we have 

d\s{t)) = -iiJ|s(t))rft-i(T2(/f-(s(t)|ff|s(t)))2|s(t))rft+-(T(if-(s(t)|ff|s(t)))|s(t))rfW^t , 

8 2 

(27b) 

and so setting t = and projecting on (s(0)| gives 

{s{0)\d\s{t))\t=o = -i{H)sdt-\a^{H-{H)sf)sdt , (27c) 

o 

with (i/"), = (s(0)|iJ"|s(0)). Substituting Eq. (|27|) into Eq. (|27|), we thus get the first 
term in the small t expansion of |(s(0)|s(t))p — 1, 

\{smsm' = l-\<^\{H-{H),f),t + 0{e) . (28a) 

This equation gives a general formula for the stochastic suppression of the quantum Zeno 
effect, independent of any assumptions about the potential. When the general form [|l5l of 
the order term coming from the standard Schrodinger evolution is included, Eq. (|28aj) 
becomes 

MOMm' = 1 - ((^ - {H)sr)s (^\aH + t'^ + 0{aH') + 0(1') ; (28b) 

in other words, the first two terms in the small t expansion are governed to leading order in 
cr by the initial state energy variance. When the potential is assumed to obey Eq. (pa|), we 
have 

{H)s = Es + Vss = Eg , 

{{H - {H%r)s = {H')s - {H)l = iV%s , (28c) 
and so Eq. (|28b|) becomes 

|(.(0)|.(t))p = 1 - {V'),, (^^ah + t'^ + 0{aH') + 0{t') , (28d) 
in agreement with the result of Eqs. (pT|a,b) and the unitarity sum rule. 
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X. DISCUSSION AND ESTIMATES OF BOUNDS ON a 



We have seen that to leading order in the perturbing potential, the stochastic terms 
governed by a do not affect either the Lorentzian line profile or the transition rate per unit 
time as evaluated in the Weisskopf-Wigner approximation, but only produce a change in the 
short time transient behavior of the transition probabilities from the initial state. This is 
a direct result of the fact that the energy-driven stochastic Schrodinger equation is energy 
conserving. On dimensional grounds, the transition rate per unit time F could contain, in 
addition to the usual terms of the form 5{Es — Em)\Vms\'^ ■, a term of the form cr^(V"^)ss- 
However, this additional term is not energy conserving, and as a result we have seen that 
its coefficient precisely cancels to zero in the Weisskopf-Wigner approximation. 

Because the transition rate per unit time and Lorentzian line shape are unaffected by 
cr, bounds on a from particle decays result only from experiments in which a metastable 
system is monitored as function of time from a known time (or vertex location) of formation. 
According to Eqs. (|28|a-d), for small times the effective transition rate per unit time is 

Vn = -^a\^Ef = \a\V')ss , (29a) 

with {AEf = {{H - {H)s)^)s the initial state energy variance. This can be interpreted 
as an early time decay rate coming from spontaneous reduction induced by the stochastic 
fluctuation terms, in agreement with the estimate F/j ~ a'^{AE)'^ used in earlier discussions 
0, 1^. In order for the rate of Eq. ( |29a| ) to not lead to pronounced early time deviations 



from the observed decay rate F, we must have 

Tr<T , (29b) 
which writing o"^ = M^^ implies the bound 

iV^)ss ^ Em^s l^sml^ ^ Ed . . 

4F 87TZm^jV,^ME^-E,) Svr ' ^ ^ 

with Ed deflning an energy characteristic of the decay process. In a particle physics context, 
a flrst guess would be to estimate Ed as being of order the mass of the decaying particle. 
The most massive decays for which F has been measured by tracking a metastable system 
from the point of formation appear to be vr*^ 77 decay, with an initial mass order 140 
MeV, and charmed meson decays, with an initial mass of around 2 GeV. Estimating in 
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Eq. ( pOD as the decaying particle mass, these give respective bounds on M^- of order 6 MeV 
and 80 MeV, respectively. If were significantly larger than these bounds, one would have 
observed anomalous accumulations of decay events close to the production vertex, as a result 
of decays induced by spontaneous reduction. For comparison, the observation of coherent 
superpositions of energy eigenstates in the neutrino, if-meson and 5-meson systems gives 
bounds ig, respectively, of M,, > lO-^OQeV, > 2 x IQ-^^GeV, and > 2 x l{]-^^GeY. 

Thus the charmed meson decay bound on represents a significant improvement over 
the coherent oscillation bounds. However, it is still smaller than the Planck mass, which is 
very likely the expected value of Mo-, by a factor of 10^°! We conclude that the theory of 
decaying states in the energy-driven stochastic Schrodinger equation places only very weak 
empirical bounds on the magnitude of the stochasticity parameter a. 

We leave for future study two issues that can be addressed within the general framework 
established here. The first is an analysis of the nature of the transition between the short- 
time regime with decay rate Tr-, and the exponential decay regime with decay rate F. This is 
governed by the solution of Eqs. (p!7|a,b) before making the Weisskopf-Wigner approximation 
of replacing E in the denominator of the final term in Kab by Es- The second is an analysis 
of the magnitude of the energy defined by Eq. (^), for various dynamical models of the 
decay process, as reflected in the energy spectrum of the unperturbed states \m) and in the 
magnitudes of the decay- inducing matrix elements Km- 
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XII. ADDED NOTE 



Lajos Diosi jl^ has pointed out an elegant stochastic-theoretic technique that allows the 



main physical results of this paper to be derived in a few lines, starting from the standard 
quantum mechanical results that hold when the stochasticity parameter a is zero. Diosi 
makes three principal observations. The first is that the quantities of direct physical interest, 
as pointed out in Sec. [V|, are the expectations £'[|Cm.(t)P] of the squared magnitudes of the 
perturbation coefficients. Since according to Eq. ( |BB| ) we have 

|C^(t)P = \{mm))\^ = {m\m){mH = {m\pit)\m) , (31a) 
with p{t) the density matrix 

pit) = m)){m\ , (31b) 

to calculate E[\Cm{t)\^ it suffices to know E[p(t)], in other words 

E[\CUt)\'] = {m\E[p{t)]\m) = Tr(|m)(m|)E[p(t)] . (31c) 

Since the dynamics of E[p(t)] is governed by the Lindblad-type equation of Eq. (|^), to 
calculate the physically relevant expectations it thus suffices to solve the dynamical problem 
specified by Eq. (||), supplemented by the initial condition 

E[p{0)]=p{0) = \sa){sa\ . (31d) 

Diosi's second observation is that the dynamical problem specified by Eq. (^, with the 
initial condition of Eq. ( |31d|) , can be compactly solved by a simple stochastic trick. The trick 
uses the fact that there is a second stochastic Schrodinger equation, simpler in structure 



than that of Eq. (p!aD, which also leads to Eq. (^) as the equation for the evolution of the 



stochastic expectation of its density matrix. To see this, consider the stochastic Schrodinger 
equation 

d\ij) = -iH\ip)dt - -a^H^\ij)dt + -iaH\il))dWt . (32a) 
8 2 

This equation differs from that of Eq. (|1^) in having an imaginary noise term, with operator 
coefficient if, instead of a real noise term with operator coefficient H — (H). A simple 
calculation, using the Ito calculus rules of Eq. (|lBD, shows that Eq. ( |32a|) also leads to 



preservation of the norm of the state {ip), and leads to the density matrix evolution equation 

dp = i[p, H]dt - la'^iH, [H, p]]dt + lia[H, p]dWt , (32b) 
8 2 
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which has the stochastic expectation 

1 



dE[p\ = i[E[pl H]dt - -a'[H, [H, E[p]]]dt , (32c) 



which is identical to Eq. (^). Hence the imaginary noise equation of Eq. ( |32a{ ) will lead to 
the same results for the physical quantities ii^[|Cm(t)p] as the real noise equation of Eq. 
even though the stochastic details of the two processes differ! 



Diosi's third observation is the fact that Eq. (p2aD can be immediately formally integrated 
to give 

\m) = exp[-iH{t - ^aWt)]\m , (33a) 



as can be readily ascertained by use of Eq. (|2q ) with the choice 

a = UaH . (33b) 

Combining this observation with the first two, then leads to a very simple rule for calculating 
the stochastic modifications of decay processes governed by Eq. ([Ta|). Let i?[|C^(t)p] be the 
quantities of physical interest, viewed as functions of a as well as oft, so that E[\C^^{t)\^] = 
|C^(t)P are their values as calculated from the standard Schrodinger evolution with no 
stochasticity. Then Eqs. ( |31a| ) through (33b) imply the simple relation 

E[\C:^m = E[\C'^{t - ^aW,)\'] , (34a) 

between the probabilities calculated in the standard Schrodinger analysis, and the stochastic 
expectations of the probabilities as calculated from Eq. (0). The recipe is simply this: take 
the known expressions for the probabilities calculated in standard quantum mechanics, re- 
place t by t — ^crWt, and take the stochastic expectation. The needed stochastic expectations 
of powers of Wt can all be read off from the expansion of Eq. (|3c|) in powers of a, 



E[Wt\ = , E[W^] = t , E[Wj^] = , =3t\... . (34b) 

Let us now apply Diosi's observations to rederive the principal results found above for 
the stochastic analog of the Weisskopf-Wigner analysis. First, let us consider the short time 
behavior of the survival probability given in Eq. ( p8b| ). The standard answer when a = 0, 
which gives the quantum Zeno effect, is 

\{smsm'=l-{{H-{H),y),t' + 0{t') . (35a) 
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Following the recipe, we have 



E[{t - ^(rWt f] = E[t^ - taWt + ^a'^W^] = + ^ah . (35b) 



On substitution into Eq. ( |35a| ) this gives for general a 



MOMm' = l-{iH- {H),r),{t' + \aH) + (35c) 

in agreement with the result for the stochastic modification of the quantum Zeno effect given 
in Eq. (pSB]) . 



Next let us apply the recipe to the formula for the initial state survival probability 
obtained using the Weisskopf-Wigner approximation, which is valid for times t that are not 
too small (and also not too large). The standard analysis gives 

|C°(t)p = exp(-n) , (36a) 

with r the Golden Rule decay rate of Eq. ( |19b| ). Replacing thy t — \crWt and using Eq. ( ^c] ) 
to take the stochastic expectation, we get as the exact formula for the stochastic modification 
of the Weisskopf-Wigner approximation 

^[|CrWr]=exp[-r(l- Vr)t] , (36b) 

which reduces, when the correction term of relative order a^F is neglected, to the answer 
found in Eq. (|20b| ) . Since cr^F ~ 0(cr^V^^), we see that the calculation of Sees. 4-6 above 
did not succeed in keeping all terms of order a'^V'^, and in fact there is a small stochastic 
correction to the decay rate, with the corrected decay rate given by 

F'^ = F(l - -a^F) . (36c) 
8 

However, writing o"^ = as in Sec. ^ as long as > Eg this correction is not significant 
within the Weisskopf-Wigner approximation scheme, which treats the line width F as a small 
quantity relative to Eg. 

Finally, let us apply the recipe to the formula giving the probability for a transition to 
the state \m). The standard Weisskopf-Wigner approximation result for this is given by 
Eq. (pOcD with a = 0, 



|2 



- 2exp[-iFt] cos[(E, -E^ + M)t] ) . (37a) 
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Applying the recipe, and again using Eq. (|3^) to evaluate the needed expectations, we get 
the exact stochastic extension of Eq. (p7a| ), 

E\\CZ,^sit)n = 7 [l^kL ^(^exp[-r(l--a2r)t] + l 

- 2exp[-ir(l - - la\E, - + Mft] 

z io o 

X cos[(E,-E„ + M)(l- Vr)t]^ . (37b) 

Again, when simplified to leading order in V , this gives the result of Eq. ( pOc| ) above. 
However, even before dropping nonleading terms in V , we see that Eq. ( |37b|) implies the 
Lorentzian formula of Eq. ( P0d|) in the large time limit. 

From the above exposition, we see that Diosi's observations not only greatly simplify the 
calculation of the physically relevant quantities, but also give results that are completely 
independent of the assumption of Eq. (pa]) that was used to linearize the stochastic equation. 
(This is something that one might have already suspected from the fact that Eq. ( p8b| ) is more 
general than Eq. ( |28d| ).) Thus, the only approximations that are needed to get stochastic 
results are those that are used in the standard, non-stochastic quantum mechanical analysis. 
Moreover, the "miraculous" cancellation of the cx^ terms in the Weisskopf-Wigner approx- 
imation to the mass and decay matrices, exhibited above in Eq. (|18a]) , is given a deeper 
explanation. There is an extensive literature |jl9| discussing the decay problem without 
making the Weisskopf-Wigner approximation (i.e., without replacing E by E^ in the order 
V'^ terms of the Laplace transform kernel Kah)-, and these discussions can all be converted to 
results for ii^[|C^(t)p] in the stochastic case, by using the recipe of replacing thy t — \(jWt 
in the corresponding formula for |C^(t)p and taking a stochastic average over Wt. 

The relation of Eq. ( |34a| ) between stochastic and standard quantum mechanical prob- 
abilities can be applied to other problems as well. For example, the density matrix of a 
two-level system can be represented in the form 

p=^(l-^-f), (38a) 

with r = (Ti,r2,r3) the standard Pauli matrices, and with R = {Ri, R2, R3) a vector sum- 
marizing the structure of the traceless part of the density matrix. The standard, cr = 
Schrodinger equation describing Rabi oscillations of the two-level system under the influ- 
ence of an applied field oscillating at the frequency of the level separation (in co-rotating 
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coordinates, neglecting the counter-rotating field component) gives for the equation of mo- 
tion ||2^ of the vector R\a=o = 

^ = uxR\ (38b) 
at 

with = the angular frequency of precession of R^. The probabilities for finding the 
system in the upper and lower levels are given, as a function of time, by 

Pi{t) = 111 ±Rm . (38c) 



Since the general solution of Eq. (|38bD has the form 



R°{t) = Vi cos Qt + V2 sin Qt , (39a) 

with Vi 2 fixed vectors that depend on the initial state and the structure of the Hamiltonian, 
and since 

ElcosQ{t crp^f)] = exp( — cos ilt , 

2 8 

£'[sinfi(t crPVf)] = exp ( — n^(TH)smnt, (39b) 

2 8 

we have under the stochastic evolution of Eq. (p^) 

ElR^it)] = exp{--n^crH)R\t) . (39c) 



By Eq. (|38c| ), this gives for the expected probabilities when the system evolves under the 



stochastic Schrodinger equation. 



1 - ElPm = exp(-lfiV^t)[i - P°(t)] . (39d) 



This can be applied, for example, to the quantum Zeno effect experiment of Itano et. al. pi 



who carry out a proposal of Cook |22| to make repeated measurements of a two- level system 
while the vector R is precessing for a time interval t = tt/Q, for which the exponential 
damping factor in Eq. (|39c| ) becomes exp(— Ivrfio"^). Corresponding to the experimental 
value Q = 320.7 MHz and the fact that probabilities were observed to an accuracy of about 
.02 in this experiment, and were found to agree with the standard Schrodinger theory, we 
get a bound on M„ = of > 2 x 10^^^ GeV, comparable to that obtained from 

oscillations in the i^T-meson system. 
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APPENDIX 



We give here the details of the integration of the stochastic differential equation that 
appears in Sec. |V111|. Consider the linear stochastic differential equation 



dCt = {AdWt + Btdt)Ct + PtdWt + Qtdt , (A.l) 

which is to be solved for the unknown stochastic function Ct given the known functions At, 
Bt, Pt, and Qt- Although we shall proceed as if these known functions were deterministic, 
in fact all our manipulations and the final solution are unchanged p3| if the input functions 
are themselves stochastic. To solve Eq. ( |A.1|) , we transpose the Ct term on the right to the 
left and multiply by a factor Ft, which is to be determined, giving 

Ft[dCt - {AtdWt + Btdt)Ct] = Ft[PtdWt + Qtdt] . (A.2) 



We now look for an Ft which makes the left hand side of Eq. ( [A. 21 ) a total differential, up 



to terms independent of Ct that are of the same form as the terms on the right hand side. 
Making the Ansatz 

Ft = exp [ / {audWu + Pudu)] , (A.3) 
Jo 

we find by use of Eqs. and (|2^) of the text that 

d{FtCt) = Ft[dCt + atdWtCt + (A + \c^,)dtCt + atdWtdCt] , (A.4) 

which on substituting Eq. ( |A.1| ) for the final dCt on the right, and using Eq. ([T^) of the 
text, gives 

d{FtCt) = Ft[dCt + atdWtCt + (A + + (ytAt)dtCt + atPtdt] . (A.5a) 



Hence if we choose 



at = -At, l3t = -Bt + ]^Al , (A.5b) 



then Eq. ( |A.5a| ) takes the form 

d{FtCt) = Ft[dCt-{AtdWt + Btdt)Ct- AtPtdt] , (A.6) 

which by use of Eq. ( |A.2| ) becomes 

d{FtCt) = Ft[PtdWt + {Qt - AtPt)dt] . (A.7) 
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The dependence on the unknown function Ct is now entirely in the form of an exact differ- 
ential, and so Eq. ( |A.7| ) can be immediately integrated to give 

Ct = F-^ [Co + / duF^{P^dWu + {Qu - A^P^)du)] 
Jo 

/•* 1 

= exp [ / {A^dWu + (5„ - -Al)du)] (A.8) 
Jo ^ 

X (c'o + ^ c^wexp [ - ^ [A^dW^ + [B^ - ]^Al)dv)\ {P^dWu + (Q„ - A^P^)du\\j , 

which is the general solution of Eq. ( |A.1| ). 

In Sec. |V111| , we need only the case of Eq. ( A.l ) in which At = A, Bt = B, Pt = Pft, 
and Qt = Qft, with A, B, P, Q constants and with ft of the form ft = exp{Kt), and so the 
solution of Eq. ( |A.8| ) then becomes 

Ct = exp {AWt + {B- ]^A^)t) 

X + du exp ( - AWu + {K - B + ^A^)u) [PdWu + {Q - AP)du]^ . (A.9) 
Using the identity (proved by the same methods used to find the integrating factor Ft), 



exp [aWu +{13 + K)u) [PdWu + (Q + aP)d' 



u 



(A.IO) 



1 



P P O JL 

— rfexp iaWu + W + K)u) {(3 + K-i-a a^) exp iaWu + W + K)u)du, 

a a P 2 



and taking a = —A and (3 = —B + A'^/2, the PdWu term in Eq. (|A.9|) can be eliminated. 
This gives an alternate form for the solution Ct, 

Ct = exp {AWt + (5 - ^A^)t) 

X (^Co~j[exp{-AWt + iK-B + ^A^)t)-l] 

^{K-B) + Q [ duexp{- AWu + {K - B + -A^)u)^ . 
A } Jo 2 / 



+ 



(A.ll) 



Taking 



A 



E ] 



P = \(TVm.. , 



Q i^msfm 1 fm 1 „ (^E^yi Eg 

o 

K = t{Em -Es-M)-^T , 



(A.12) 
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in Eqs. ( |A.9| ) and ( |A.11| ) gives the results quoted respectively in Eqs. ( P5b|) and ( P5c|) of the 
text. 
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